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Abstract 

The thermodynamics of the infinite-range Ising spin glass with p-spin 
interactions in the presence of an external magnetic field h is investigated 
analytically using the replica method. We give emphasis to the analysis of the 
transition between the replica symmetric and the one-step replica symmetry 
breaking regimes. In particular, we derive analytical conditions for the onset 
of the continuous transition, as well as for the location of the tricritical point 
at which the transition between those two regimes becomes discontinuous. 

Short Title: p-spin Ising model 
PACS: 87.10+e, 64.60Cn 



1 Introduction 



Although the thermodynamics of the Ising spin glass with infinite-range inter- 
actions, so-called Sherrington-Kirkpatrick (SK) model Q, has been thoroughly 
investigated in the last two decades (2|, [| , comparatively little attention has been 
given to the analysis of a natural generalization of the SK model, namely, the 
p-spin interactions Ising spin glass. This model is described by the Hamiltonian 

31 



Wp(S) 



E 

l<n<i2---<i P <A r 



J, 



ip Si L Si 2 . . . S{ p 



(1) 



where Si = = 1, . . . , N are Ising spins and h is the external magnetic field. 
Here the coupling strengths are statistically independent random variables with a 
Gaussian distribution 
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Besides the acknowledged importance of the p-spin interactions Ising spin glass 
in the framework of the traditional statistical mechanics of disordered systems (it 
yields the celebrated random energy model in the limit p — > oo and the SK 
model for p = 2), it also plays a significant role in the study of adaptive walks in 
rugged fitness landscapes within the research program championed by Kauffman 

The thermodynamics of the SK model {p = 2) as well as that of the random 
energy model (p — > oo) are now well understood. In particular, for p = 2 and h = 
the order parameter function q(x) tends to zero continuously as the temperature 

(2) 

approaches the critical value T t ■ = 1 at which the transition between the spin- 
glass and the high temperature (disordered) phases takes place [||, ||. For p — > oo 
and h = 0, the system has a critical temperature T c 



at which 

it freezes completely into the ground state: q(x) is a step function with values 
zero and one, and with a break point at x = T/Tc 00 ^ ||, ||]. These results are not 
affected qualitatively by the presence of a non-zero magnetic field. In particular, 
for the SK model the critical temperature decreases monotonically with increasing 
h while the transition remains continuous, in the sense that q(x) is continuous at 
the transition line |9|, In contrast, for the random energy model the critical 
temperature increases with increasing h while the discontinuity in the step function 
q(x) decreases with increasing h and vanishes in the limit h — > oo [||, |j. 

The situation for finite p > 2 is considerably more complicated and so the 
thermodynamics of the p-spin model has been investigated for h = only [|ll], 12|. 
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In this case there is a transition from the disordered phase to a partially frozen 
phase characterized by a step function q{x) with values zero and q\ < 1. As 
the temperature is lowered further, a second transition occurs, leading to a phase 
described by a continuous order parameter function [11, 12]. However, there are 
evidences that the presence of a non-zero magnetic field decreases the size of the 
discontinuity of the order parameter q(x) leading, eventually, to a continuous phase 
transition. In fact, a recent analysis of the typical overlap q between pairs of 
metastable states with energy density e indicates that q is a discontinuous function 
of e for p > 2, and that the size of the jump in q increases with p and decreases 
with h, vanishing at finite values of the magnetic field |13{| . Moreover, a similar 
effect has already been observed in the thermodynamics analysis of the spherical 
p-spin interaction spin-glass model [|l4j]. It is interesting to note that the spin-glass 
phase of this continuous spin model is described exactly by a step order parameter 
function, i.e., the one-step replica symmetry breaking is the most general solution 
within the Parisi scheme of replica symmetry breaking 14]. 

In this paper we use the replica method to study the thermodynamics of the 
Ising p-spin interaction spin-glass model in the presence of the magnetic field h. We 
focus on the effects of h on the transition between the replica symmetric (RS) and 
the one-step replica symmetry breaking (1 RSB) regimes. In particular, we show 



that for p > 2 the discontinuous transition reported in previous analyses [11, 12 [ 



turns into a continuous one for h larger than a certain value Kt- Moreover, we 
derive analytical conditions to determine the location of the continuous transi- 
tion line, as well as that of the tricritical point at which the transition becomes 
discontinuous. 

The remainder of the paper is organized as follows. In Sec. [2] we discuss the 
replica formulation and present the formal equation for the average free-energy 
density, which is then rewritten using the RS and the 1 RSB ansdtze. These results 
are discussed very briefly since their derivations are given in detail in Gardner's 
paper [plj] . We also present the solution of the 1 RSB saddle-point equations in 
the limit of large p, extending thus the series expansions results for non-zero h. In 
Sec. H we derive analytical conditions for locating the continuous transition and 
the tricritical point between the RS and 1 RSB regimes, and present the phase 
diagrams in the plane (T, h) . Finally, some concluding remarks are presented in 
Sec. |. 
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2 The replica formulation 



We are interested in the evaluation of the average free-energy density / defined by 
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where 

Z = Tr s exp [-/3H P (S)] , (4) 

is the partition function and (5 is the inverse temperature. Here ((. . .}} stands for 
the average over the coupling strengths, and Trg denotes the summation over the 
2^ states of the system. As usual, the evaluation of the quenched average in Eq. 
(||) can be effectuated through the replica method: using the identity 



<0nZ)) = lim± In 
n— >o n 



(5) 



we first calculate {{Z n )) for integer n, i.e. Z n = n"=l anc ^ t nen analytically 
continue to n = [|- The final result is simply [11] 
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The extremum in Eq. (o) is taken over the physical order parameter 



N 



Qab = {{^J2(^)T(^)T)) a<b, 



1=1 



(8) 



which measures the overlap between two different equilibrium states S a and S b , 
and over its corresponding Lagrange multiplier \ a b. Here, (■ ■ -}t stands for a 
thermal average. To proceed further, next we consider two standard ansdtze for 
the structure of the saddle-point parameters. 
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2.1 Replica symmetric solution 



In this case we assume that the saddle-point parameters are symmetric under 
permutations of the replica indices, i.e., q a b = q and X a b = A. With this prescription 
the evaluation of Eq. (||) is straightforward, resulting in the replica-symmetric free- 
energy density 

1 1 r°° 

-Pfrs = -^ 2 \(l-q) + -(3 2 (l-q p )+ £>zln2cosh[/3E,(z)] (9) 

E s = zV\ + h (10) 



where 
and 



is the Gaussian measure. The saddle-point equations df rs /dq = and df rs /d\ = 
yield 



_P nP -i 

and 



A = Zf-i (12) 
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respectively. The replica-symmetric solution is locally stable wherever the Almeida- 
Thouless condition [0], which in this case is given by |TT|| 

A t°° 

1 - P 2 (p - 1) - / Dz sech 4 \3E S (z)] > 0, (14) 
q J-oo 

is satisfied. In fact, since Eq. (|l3| ) has either one or three positive solutions, this 
stability condition is very useful to single out the physical one. In particular, for 
h = the only stable solution is q = 0. 



2.2 Replica symmetry broken solution 

Following Parisi's scheme Q, we carry out the first step of replica symmetry break- 
ing by dividing the n replicas into n/m groups of m replicas and setting q a f, = q\, 
A a fc = Ai if a and b belong to the same group and q a b = qo, X a t = Ao otherwise. 
The physical meaning of the saddle-point parameters is the following 

1 N 

90 = <<T7E<^M^M> a<b, (15) 
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<& = «^5>?>r» , (16) 

and m = 1 — X) a • H ence 90 is the overlap between a pair of different equilibrium 
states, qi is the overlap of an equilibrium state with itself (gi > qo), and m is the 
probability of finding two copies of the system in two different states (P a is just 
the Gibbs probability measure for the state S a ). We note that in the limit n — > 0, 
the parameter m is constrained to the range < m < 1. Using this prescription, 
Eq. @ becomes 

- Pfrsb = ~^(3 2 [2Ai (1 - q\ + mq\) - 2mq Xo - 1 + (1 - m) q\ + mcft] 

i poo poo 

+ ln2 + — / Dz ln D Zl cosh m (3E (17) 

771 J —oo J — oo 

where 



H = zi vAi - A + zoVAo + h. (18) 
The saddle-point equations df rs b/dqk = yield 

h = \lV (19) 
for k = 0, 1. The saddle-point parameters q^i <?i an d m are given by the equations 

/oo 
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where we have introduced the notation 

n^xG..) cosh™/?S 

J^Dzicosh™/?- ' 1 j 

It is clear from these equations that the replica-symmetric saddle-point qo = qi = q 
is a solution for any value of m. In general, however, the 1 RSB equations will 
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admit a different solution. In particular, in the limit p — > oo the solution is 
qo = tanh 2 ((3mh), q± = I and m = c /(3 where (3 C = l/Tc°°^ is the solution of the 
equation 

(24) 



\j3 2 = In 2 cosh j3 c h — (3 c h tanh j3 c h. 



Below Tc°°\ the entropy vanishes and m sticks to its maximum value, namely 
m = 1, signaling the existence of a frozen phase in accordance with the physical 
meaning of m mentioned before. It is instructive to consider the finite p corrections 
to the infinite-p solution by expanding the 1 RSB equations around that solution. 
Thus, extending the results of Gardner [11] for non-zero h, we find 



qo = tanh 2 ((3mh) 



1 + 2£ m sech {(5mh) 



-f} 2 m 2 p/4: 
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(27) 
(28) 

(29) 



otherwise. Here, 



Y oo 



m 



1 



2i — m 



--L /°° [2 cosh (mz) - 2 m cosh m (z) 
\/2tt J—oo 



(30) 



where we have used the extended definition of the binomial coefficient to real m 



II- 

At this stage we already can realize the existence of two solutions of a quite 
different nature, signaling then the non-trivial role played by the magnetic field in 
the thermodynamics of the p-spin model. 
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A quite interesting property of the 1 RSB solution, which can easily be verified 
numerically, is that qo = for h = and p > 2, indicating thus that the equilibrium 
states are completely uncorrelated. Moreover, this result has greatly facilitated 
both the numerical and analytical analyses of the model, since the integrals over 
zq in Eqs. (p0|)-(p2]) can be can carried out trivially in that case [JO], [i2| . However, 
as explicitly shown by Eq. (p5|) the non-zero magnetic field induces correlations 
between different equilibrium states so that qo is no longer zero in this case. 

For p = 3, we present in Figs. 1, 2 and 3 the temperature dependence of the 
RS and 1 RSB saddle-point parameters for h = 0, 0.5 and 1, respectively. As 

mentioned before, for ft = we find q = qo = 0. The size of the jump in qi 

(3) 

decreases with increasing h and disappears altogether for h > K* T ' 0.57. Of 
particular interest is the temperature dependence of the saddle-point parameter 
m: at the discontinuous transition it reaches its maximal value, namely, m = 1, 
while at the continuous transition it assumes a certain value m = m c < 1, which 
depends on T and p. As expected, the behavior pattern depicted in Fig. 3 is very 
similar to that found in the analysis of the magnetic properties of the SK model 
flC)[| , as the transition is continuous in that model. We note that since m plays no 
role in the RS solution, the curve for m must end at the transition lines. 

The location of the transition lines as well as the characterization of the critical 
values of the saddle-point parameters are discussed in detail in the next section. 



3 Transition lines 

As indicated in the figures presented before, there are two qualitatively different 
types of transition between the RS and the 1 RSB regimes which we will discuss 
separately in the sequel. 

3.1 Continuous transition line 

The location of the continuous transition between the RS and the 1 RSB solution 
is determined by solving the 1 RSB equations in the limit of small q\ — qo- More 
pointedly, subtracting Eq. (pC] ) from Eq. (^|) and keeping terms up to the order 
(<7i - 90 ) 2 yields 

qi ~ q °- PHP-1)X B^n) (31) 

where 

\ roo 

Bo (qo) = 1 - P 2 (P ~ 1) — / Dz sech 4 [(3E (z)] , (32) 

Q0 J-oo 
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and 



B 2 (qo,m) 



p-2 + 4(3 2 (p-l) A (3 - 2m) 



Dz sech 4 [pEo (z)] 
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(33) 



Here 



H (z) = zy/Xo + h (34) 

with Ao given by Eq. (fjl|). We note that both Bq and B 2 are negative quantities in 
the 1 RSB regime. Since at the continuous transition qi — > qo —* where q is the 
replica symmetric saddle-point (|l^), the transition line is given by the condition 



Bo (q) = 



(35) 



which, as expected, coincides with the replica-symmetric stability line given by Eq. 
(|i"4|). To specify the value of m at the critical line, denoted by m c , we expand Eq. 
(p2]) for small q\ — qo (in this case we must keep terms up to the order (gi — qo) 3 ) 
and then subtract it from Eq. (20). Using the condition (|35| ) together with qo ^ q 
yields 

B 2 (q, 1) 



B±{q) 



(36) 



where 



/oo 
Dz sech 2 [/3Z S (z)} tanh 2 [/3H S (z)} 
-oo 

/oo 
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-oo 



(37) 



with E s and A given by Eqs. ([h]) and (12), respectively. Eq. (|36|) holds provided 
that m c < 1 and so the continuous transition line must end at a tricritical point, 
whose location is obtained by solving 



B 2 (q, 1) = 



(38) 



and Eq. (^) simultaneously. As usual, the denominator in Eq. (^) vanishes at 
the tricritical point. 



3.2 Discontinuous transition line 



The location of the discontinuous transition line is determined by equating the free 
energies of the RS and 1 RSB solutions, given by Eqs. (g) and (17), respectively. 
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This task is greatly facilitated in this case by noting that setting m = 1 in Eq. (pC]) 
yields qo = q. Moreover, since for m = 1 Eq. ([l7]) becomes independent of q\ (and 
Ai) one has frsbim = 1) = frs- Thus, for fixed h the temperature at which the 
discontinuous transition takes place is obtained by solving the 1 RSB saddle-point 
equations with m = 1 for q\, qo = q, and T = T c . 

3.3 Analysis of the results 

The phase diagrams in the plane (T, h) are presented in Figs. 4, 5, and 6 for 
p = 2, 3 and 10, respectively. The solid curves are the RS stability condition, Eq. 
(|j~i|), whose upper branch coincides with the continuous transition line, Eq. (|35|), 
The discontinuous transition lines (dot-dashed curves) join the continuous ones at 
the tricritical points (full circles). We also present the lines at which the entropy 
of the RS solution vanishes (short-dashed curves), which for h = intersect the 
temperature axis at T = 1/ ^2\/\n 2j ~ 0.60, whatever the value of p > 2. We 
note that for p — > oo the condition for the vanishing of the RS entropy yields 
exactly the discontinuous transition line for the random energy model, Eq. (pi|). 
The agreement between these lines is already very good for p = 10 and h not too 
near ht, as illustrated in Fig. 6. 

Since our results are valid for real, though physically meaningless, values of 
p > 2 as well, we present in Fig. 7 the value of the saddle-point parameter m 
at the continuous transition line, given by Eq. (|36|), for several values of p. As 
expected, for p > 2 we find m c = 1 at the tricritical points. We mention that, 
despite the numerous studies of the SK model, we are not aware of any calculation 
of the 1 RSB saddle-point parameter m over the Almeida-Thouless stability line. In 
Figs. 8 and 9 we present the values of T and h at the tricritical point, respectively, 
as functions of the real variable p. For p — » 2 we find ht — > and Tj — ► 1, while 
for large p we find that ht increases like y/p\np, Tt like \Jpj lnp and 1 — qt goes to 
zero like 1 / (p^/ln p) . These results indicate that the phase diagrams in the plane 
(T, h) do not display the two types of transitions only in the extreme cases p = 2 
and p — ► oo. 



4 Conclusion 

Some comments regarding the validity of the 1 RSB solution are in order. The 
stability analysis of that solution carried out for h = indicates that it becomes 
unstable for low temperatures [11]. A seemly simpler approach to check the phys- 
ical soundness of the 1 RSB solution is to evaluate numerically its entropy. This 
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procedure, however, has proved very elusive: since the entropy becomes negative 
when it is of order e~^ p , the numerical precision required to determine the tem- 
perature at which it vanishes T" is exceedingly large. For instance, for h = we 
find T" = 0.10, 0.087 (0.19) and 0.034 (0.18) for p = 2, 3 and 5, respectively. 
The numbers between parentheses are the numerical estimates of ref. |i~2|1 . As our 
numerical results are in good agreement with that of ref. |l6[ for p = 2, and also 
are consistent with the trend of decreasing T" with increasing p, we think they 
are the correct ones. Already for p > 5, however, we have failed to obtain reliable 
estimates for T" . Since the precision problem becomes much worse for non-zero h, 
due to the numerical evaluation of the double integrals, we refrain from presenting 
the estimates for T" in that case. 

Although for finite p the 1 RSB solution certainly does not describe correctly 
the low temperature phase of the p-spin Ising spin glass, it probably yields the 
correct solution near the transition line delimiting the replica symmetric and the 



replica symmetry breaking regimes. In fact, according to Gardner 11], considering 
further steps of replica symmetry breaking within Parisi's scheme will result in a 
new continuous transition between the 1 RSB regime and a more complex regime, 
described by a continuous order parameter function. In this sense, we think that 
our results regarding the transition lines between the RS and the 1 RSB regimes 
are not mere artifacts of the replica method but indeed describe genuine features 
of the thermodynamics of the infinite range p-spin Ising spin glass in a magnetic 
field. 
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Figure captions 



Fig. 1 One-step replica symmetry breaking saddle-point parameters m (solid 
curve) and q\ (short-dashed curve) as a function of the temperature T for p = 3 
and h = 0. In this case qo = q = 0. The discontinuous transition occurs at 
T « 0.65. 

Fig. 2 One-step replica symmetry breaking saddle-point parameters m (solid 
curve), qo (long-dashed curve), q\ (short-dashed curve) as a function of the tem- 
perature T for p = 3 and /i = 0.5. The dash-dotted curve is the replica symmetric 
saddle-point parameter q. The discontinuous transition occurs at T « 0.74. 

Fig. 3 Same as Fig. 2 but for h = 1. The continuous transition occurs at T « 0.66, 
at which m « 0.87. 

Fig. 4 Phase diagram in the plane (T, h) for p = 2. The RS saddle-point is 
locally unstable inside the region delimited by the solid curve, which coincides 
with the continuous transition line between the RS and the 1 RSB regimes. The 
short-dashed curve delimits the region inside which the RS entropy is negative. 

Fig. 5 Same as Fig. 4 but for p = 3. The RS stability line coincides with the 
continuous transition line in the branch above the tricritical point (full circle), 
located at T t rs 0.74 and h t ~ 0.57. The dot-dashed curve is the discontinuous 
transition line. The convention is the same as for Fig. 4. 

Fig. 6 Same as Fig. 5 but for p = 10. The tricritical point (full circle) is located 
at T t w 1.01 and h t « 3.07. 

Fig. 7 Saddle-point parameter m at the continuous transition line for (from 
bottom to top at T = 0.5) p = 2, 2.01, 2.1, 10 and 3. 

Fig. 8 Temperature at the tricritical point Tt as a function of p > 2. Only integer 
values of p have physical meaning. 

Fig. 9 Magnetic field at the tricritical point ht as a function of p > 2. Only 
integer values of p have physical meaning. 
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